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I. INTRODUCTION 

In a recent paper Carbotte et a/.0 showed that the 41 meV spin resonance observed in 
spin polarized neutron scattering experimentJi on YBa2Cu306.95 (YBCO) has a counterpart 
in the infrared optical conductivity in which an optical resonance is seen at the same energy 
demonstrating its coupling to the charge carriers.iS The temperature evolution of the spin 
resonance can be determined from consideration of the frequency dependence of the optical 
scattering rate measured at several different temperaturedl in the superconducting state. 
The required analysis which allows one to extract the resonance from the optics is carried 
out within a generalized Eliashberg formalism for a d-wave superconductor.! The end re- 
sult is a charge carrier-exchange boson excitation spectral density which characterizes the 
inelastic scattering. The spectral density /^x(co') shows a distinct evolution with decreasing 
temperature. At Tc the spectrum obtained is linear in u at small u, followed by broad peak 
at some characteristic energy ujsf and then there is a very slow decay at higher u extending 
up to a cutoff ujc of order 400 meV. This spectrum can be fit with the spin fluctuation form 
employed in Pines' groupiil i+(i^^sp)i ^ i'^c — ^) (MMP model), where P is a coupling 
constant to the charge carriers. As the temperature is lowered below Tc, the fluctuation 
spectrum obtained from the infrared data undergoes important changes. It develops a gap 
at the lowest energies as well as a peak at 41 meV. The strength of this peak tracks well the 
observed growth of the area under the spin susceptibility, 53m{x(q, cj)} at q = (tt, tt), mea- 
sured by Dai et a/.i in neutron scattering. The higher energy part, however, remains largely 
unaffected. The observed changes in spectral density are as expected and are interpreted as 
coming from feedback effects on the excitation spectrum due to the changes brought about in 
the electronic system by the onset of superconductivity. This arises in any purely electronic 
mechanism of superconductivity in which the pairing proceeds through bosons exchanged 
between the charge carriers and where the bosons themselves are intrinsic to the electronic 
system.1'0 The changes in the low energy part of /^x(co') lead directly to the phenomenon 
referred to as the collapse of the quasiparticle scattering rate and is responsible for the freez- 
ing out of the inelastic scattering at low temperatures. In turn, this leads to a prominent 
peak around 30 K in the in-plane microwave conductivityS as a function of temperature T 



and a corresponding peak in the electronic thermal conductivity.lll 

Schachinger et a/.i have recently found that this same spectral density Pxi^) ^-Iso gives, 
within an Eliashberg formalism, good agreement with observed properties of the supercon- 
ducting state such as the ratio of the gap to the critical temperature, the fractional optical 
spectral weight involved in the condensation (superfiuid stiffness), the temperature depen- 
dence of the penetration depth, the magnitude of the condensation energy, as well as other 
quantities. It is also obvious from the way the temperature dependent spectral density has 
been derived in Ref. 5 that the coupling to the 41 meV spin resonance in YBCO cannot play 
the role of the 'glue' leading to superconductivitylii because it is absent at Tc. 

In this paper we consider the temperature dependence of the microwave conductivity for 
five different frequencies between 1 and 75 GHz observed recently0 in ultra pure samples 
of YBa2Cu306.99 grown in BaZnOa crucibles. Our calculations are based on the previously 
determined spectral density Px{^) for twinned YBCO single crystals which is not modified 
in any way. It is found, however, that to understand the low temperature data (10 — 20 K) 
it is necessary to additionally introduce some elastic impurity scattering. A model impurity 
potential intermediate between weak (Born) and strong (unitary) scattering is developed 
which provides a reasonable, if not perfect overall fit to the data. 

In Sec. II we give the necessary formalism. This is followed by the presentation of our 
results in Sec. III. Comparison with the data is also presented in this section. Section IV 
contains discussion and a conclusion. 



II. FORMALISM 



We consider a ci-wave superconductor. To include the inelastic scattering which is known 
to be strong in the cuprates we need to go beyond a simple BCS formalism. The minimum set 
of equations that allow us to do this are generalized Eliashberg equations. These equations 
involve two channels. The ordinary renormalization channel which remains in the normal 
state and leads to renormalization of the Matsubara frequencies by the interactions. The 
second is the pairing channel which we assume to have d-wave character. For simplicity we 
take the pairing potential to be separable in incoming (k) and outgoing (k') momentum in 
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the two dimensional Cu02 Brillouin zone. The charge carrier-fluctuation spectrum spectral 
density Pxi^) which characterizes the inelastic scattering in the renormalization channel is 
assumed to also cause the pairing and while, in principle, the spectral density could have 
a different shape in this case, for simplicity, we ignore such complications. We do allow, 
however, for a possible difference in magnitude through a constant g. For the renormalized 
frequencies a) we use Px{^) a)(a; + i5) is isotropic while for the pairing energy A we 
use gPxi^) cos(20) cos(2(/)') which immediately leads to a d-wave form A{uj + iS) ~ cos(2(/)), 
where is the angle defining the direction of momentum k on a cylindrical Fermi surface. 
The generalized Eliashberg equations which play the central role in this study are 

m=o \ ya)2(ia;„) + A2(ia;^; 0') 

+mg J dz cos{2<f>)Px{z) [n{z) + f{z - u)] 

— oo 

^ / A(z/-z + ^^;0Ocos(20O y ^^^^ 
and in the renormalization channel 



UJ 



[u + iS) — u + mT — iUm) — A(z/ + iUm)] 

oo 

+iTT / dzl\{z)[n{z)+f{z-iy)] 



I 



\J uj'^[i' — z -\- lb) — — z -\- id](j)') I ^ ' ^ ' 

In the above A(ia;m; 0) is the pairing energy evaluated at the fermionic Matsubara fre- 
quencies ujm = 7iT{2m — l),m = 0, ±1,±2, .... f{z) and n{z) are the Fermi and Bose 
distributions respectively. The renormalized Matsubara frequencies are ij{iu!m)- The ana- 
lytic continuation to real frequency i/ of the above is A(i/ -|- iS; (p) and u;{u! + iS), where 5 is 
a positive infinitesimal. The bracket (• • • ) is the angular average over 0, and 

oo 
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D(^)-{ , + ). (3) 



nM = ( . >. (4) 



Equations (|T]) are a set of nonlinear coupled equations for the renormalized pairing potential 
A(z/ + i6; (f)) and the renormalized frequencies uj{v + i5) with the gap 

A(zy + ^5;0) = ^_^-^, (5) 

where the renormalization function Z{v) was introduced in the usual way as Cj{v + i5) = 
vZ{v). In Eq. sets the size of the impurity scattering and c is a parameter related 

to the impurity potential. The Born or weak scattering limit corresponds to a large value 
of c while the unitary or strong scattering limit corresponds to c = 0. The impurity term is 
obtained from a T-matrix approach to the impurity problem.lll It does not include all possible 
complications that have come to be known as possibly of some importance in the cuprates. 
Recent STM dataSIll has revealed significant inhomogeneities and impurity studies based on 
BdG equations have shown that the superconducting order parameter is strongly modified in 
the vicinity of an impurity.lli0illii'§i@iiilil Such complications go well beyond the present 
approach. Here we will not attempt a complete microscopic description of the impurity 
scattering but instead treat r+ and c as parameters which we will determine through a 
best fit to data. To get a better fit may well require the introduction of more sophisticated 
effective potentials. 

In the pure limit, i.e.: no impurity contribution (r+ = 0) all parameters in Eqs. (0) that 
serve to characterize the particular material of interest, are fixed from our previous work. 
The charge carrier-fluctuation spectrum spectral density Pxi'^) which enters through Eq. @) 
was previously obtained by Schachinger et al% from infrared optical data using an inversion 
techniques which allows one to construct from the optical scattering rate. In the 

present calculations we simply use these results without any modifications. The parameter 
g is also fixed, and was determined to get the measured value of the critical temperature 
in YBCO. While we will first present results in the clean limit we will later see that, to 
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understand the low temperature data (T 0) we will need to consider some impurity 
scattering. The value for r+ and c will be chosen to best fit the data as we have already 
emphasized. 

The optical conductivity follows from a knowledge of uj and A. The formula to be 
evaluated is 



n 2 



X ( / dt/tanh f-^^ — ; -r ; — 

\2TJ E{u;4>) + E{u + n; 

' 



+ du tanh 



1 



2T / EHu; 











E*{iy + n;(l)] 



[1 - N{u; (p)N{v + VL](t))- P(z/; (/))P(z/ + VL; 
— [l-N\iJ-<P)N''{y + n-<P) 



tanh 



2T 
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+ / dz/tanh 



[1 + Ar^(z/; <P)N{u + VL-<P) + P\u- 0)P(z/ + Vt- , 
1 



2T J 1 EHv\(t)) + E*{v + ^ 

-n 

1 



E(u + n:d)) - E*(u:. 



1 + N*{u; (j))N{u + fi; 0) + P*{u; (j))P{u + fi; : 



with 



and 



a)k(a;) 



(6a) 



(6b) 



(6c) 



P(a;;0)' '^-'-^^ E(a;;0)" 
In the above, (■ ■ ■ ) means, as before, an average over the angle (p and the star refers to the 
complex conjugate. A^(0) is the electronic density of states at the Fermi surface and vp the 
Fermi velocity. The prefactor in (6a) can be worked out to be proportional to the plasma 
frequency squared, f^p/47r = ne^/m*. Here n is the electron density, e the electron charge, 
and m* its effective mass. 
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III. RESULTS AND DISCUSSION 

In Fig. which has five frames we show our results for the temperature dependence of 
the real part of the microwave conductivity ai{u!) in the pure limit (open triangles) i.e.: 
r+ = in Eq. (|l]) at the five measured frequencies of Hosseini et a/.0 They correspond from 
top frame to bottom frame to = 1.14, 2.25, 13.4, 22.7, and 75.3 GHz respectively. For the 
lowest frequencies considered, the agreement with the data (solid squares), which we read of 
a graph in Ref. 11, is not good but the agreement does improve as the microwave frequency 
is increases. In particular, in the two top frames the height of the calculated peak is too 
high and it falls at a somewhat lower temperature than indicated in the measured curve. 
These deficiencies can be largely removed when a small amount of impurity scattering is 
additionally included. The effect of impurity scattering will show up most prominently at 
the lowest temperature where the inelastic scattering is becoming very small. 

Before proceeding with a fit to the data which includes both impurities and the inelastic 
scattering it is useful to first consider the BCS limit of our generalized Eliashberg equations 
(0) and to understand the effect of impurities in this instance. At low temperatures, the 
inelastic scattering rate which depends on real processes is small and the impurities will 
dominate; thus the BCS theory will become more applicable although it does ignore all 
renormalizations from the inelastic interaction (virtual processes). 

To obtain the BCS equations from (|l|) we ignore the effect of Px{^) the renormalization 
channel, so 

and in the gap channel we assume that the Boson frequency in Px{^) is very high compared 
with all other energies of importance. This means that Bose and Fermi factors in the second 



term on the right hand side of ([Ta|) are negligible and we can also replace the X{i'±iuJm) by 
a constant (A) with a cutoff (c^c) on the Matsubara sum applied to get convergence. This 
gives a gap A(0) independent of frequency z/ which satisfies the equation 

A(0) = 2nTg cos(20) A / ^(fW^V . (g) 
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The impurities enter directly in Eq. and affect the gap given in Eq. (|]) through the renor- 
mahzed Matsubara frequency uj{iujjn) which appears in the square root in the denominator. 
In the hmit of z/ = we can write 

i(. + «)«n--r--J|2l^ (9) 

which can be solved self consistently for the impurity scattering rate at zero frequency. We 
can evaluate = 2K{Ao/lj)/tt where K{x) is the elliptic integral of the first kind.il The 
quantity appearing in Eq. is ^(^7) = [27/(7rAo)]ln(4Ao/7) and for a general value 

of c, the equation for 7 is 

^In ( ^) 

This equation shows that the self consistent impurity scattering rate 7 at zero frequency in 
the superconducting state is strongly dependent on the parameter c. For the strong coupling 
unitary limit c = an approximate solution has been given by Hirschfeld and Goldenfeldii 
for vrr"*" ^ Aq as 



7 ~ 0.63V7rr+Ao. (11) 

Note that 7(0 = 0) is much larger than vrr^ in this limit. In the opposite limit (Born limit 
or weak scattering potential) c 00 and ttT^/c^ is to be replaced by ttFtv and 

^(c ^ 00) = 4Aoe-^"/(2^^\ (12) 

which shows that 7(0 00) is now much smaller than the normal state value of F, Ftv- 

This analysis demonstrates that the zero frequency self consistent scattering rate in the 
superconducting state is much larger than its normal state value in the unitary limit but 
is much smaller in the Born limit. In particular, this implies that in the Born limit the 
impurity limited quasiparticle mean free path for a given impurity content will be much 
larger in the superconducting state than in the corresponding normal state if the inelastic 
scattering is ignored. 

In Fig. 1^ we illustrate the general case. What is plotted is 7(c) as a function of c for 
a specific value of F+ = 0.15 meV and a zero temperature gap of Aq = 24\/2meV. The 
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underlying normal state scattering rate with which 7(c) is to be compared is 7rr"'"/(l + c^) 
for any value of c. In this example for c = 0, 7(c) is larger than 7rr+ (by a factor of 5), 
for c = 0.2 they are comparable and for c = 0.3 it is already much less. By changing c we 
can change the value of quasiparticle scattering rate at = in the superconducting state 
by orders of magnitude and this will be of importance for our analysis of the experimental 
data. 

It is instructive to look as well at the frequency dependence of the underlying quasiparticle 
scattering rate, or more precisely, the imaginary part of the renormalized frequency, namely 

.-'M^9m.M^i.M^.r^-J||l-j. (13) 

In Fig. ^ we show results for several values of c. The behavior of t^^(z/) vs. v at small v 
changes radically with choice of c. In the unitary limit there is a small region where T^^iv) 
is fairly flat but for finite c, r^^iv) begins to look like a d-wave quasiparticle density of states 
and the scattering is radically affected by the onset of superconductivity which modifies the 
density of final states available for scattering. These effects can be understood simply in 
the clean limit F"*" and for temperatures T > 7. This limit is considered in the work of 
Hirschfeld tt a/.@ who treat the two cases c = and c 00 explicitely. Here we consider 
finite c and v small (^ Ao)0 

Aoc4 + 2cM_(z/) + (4z//Ao)2A+(z/) ^ ^ 



with 



2 



(15) 



,2/ V uj 

It is clear that for c — >■ 00 r~^(z/) becomes proportional to u while for c = it goes like u'^ 
at small u. For a general c, the quasiparticle scattering rate r~^(i/) is importantly dependent 
on u and, therefore, is quite different for the constant of the familiar normal state Drude 
model. This means that, while we have two parameters F"*" and c to adjust, the underlying 
complicated variation of r~^(z/) with u gets reflected directly in the frequency variation of 
the conductivity and leads to a non-Drude form in sharp contrast to the underlying normal 
state. 
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We present results in Fig. ^. What is plotted is the real part of the conductivity (Ti{uj, T) 
in units of A^(0)f |. as a function of frequency uj in the range up to 1.0 meV (~ 242 Ghz). The 
temperature is set at T = lOK and the impurity scattering at = 0.15 meV in Eq. (|^. 
Various values of the impurity parameter c are shown. The solid gray curve in Fig. ^ is 
the normal state shown for comparison. It displays the classical Drude form with a Drude 
width of 7rF+. The other curves are in the superconducting state at T = 10 K with the zero 
temperature (i-wave gap taken to be Aq = 24\/2 meV. The black solid curve is the Born limit, 
dotted the unitary limit, dashed c = 0.4 and dash-dotted c = 1.0. None of the curves in the 
superconducting state follow the Drude form of the normal state aiiuj) = 2Timp/[l + {ujTimp)'^] 
with Timp = 2.12 meV^^. The curves near the Born limit show a concave up rather than the 
concave down behavior of the normal state Drude. The curve for the unitary limit is much 
flatter than the solid gray curve reflecting a value of 7(c) which is large compared with its 
normal state counterpart. The inset in the top right hand shows this on a different scale and 
allows the reader to better see the radical difference in behavior between Born and unitary 
limit. Although, neither of these limits show a Drude variation with u, we can still think of 
the half width of each curve as giving a measure of the underlying quasiparticle scattering 
rate. One then concludes that in the Born limit it is much smaller than in the normal state 
while in the unitary limit it is much larger. 

The data in the work of Hosseini et a/.0 which we use here for comparison with theory 
were flt to a Drude form and the authors concluded that the effective quasiparticle scattering 



rate was fairly temperature independent and constant. It is clear from our Eq. (14) and Fig. y 
that this behavior will never be reproduced in a BCS theory with only elastic impurity 
scattering whatever the value of c. The underlying quasiparticle scattering rates are always 
highly frequency dependent and this modulates the Drude line width as u changes. This 
also implies a temperature dependence since a change in T involves a different sampling 
of the frequency dependence of t'^{lj). Any two fluid approach would need to account 
for these features of the scattering rates of the normal quasiparticles as well as the energy 
dependence in the density of states.^! Including some impurity scattering in addition to the 
inelastic scattering in Eq. greatly improves the agreement with experiment. The half 
width of the conductivity as a function of u in our clean limit Eliashberg calculations is 
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significantly smaller than what is observed experimentally. (See Fig. ^, open symbols.) This 
does allow us to add some impurities which of course always result in an increase of the half 
width. 

To make an appropriate choice of r+ and of c we are guided by the Drude analysis of 
the data provided by Hosseini et a/.0 They find a very narrow width to their Drude of the 
order of 1/3 K. It is clear that the amount of impurities involved is very small and that 
7rr+/(l + c^) = 7Ar(c) which gives the scattering in the normal state is correspondingly 
small. It is also clear that the unitary limit is unlikely since this increases 7(c) considerably 
as compared with the normal state equivalent 7iv(c) which would then have to be much 
smaller than 1/3 K which is not likely. On the other hand, the Born limit gives a concave 
upward curve which drops too rapidly as u increases out of zero. It appears that some 
intermediate c-value is favored but on examination of Fig. |^ and the form ([l^ it is clear 
that our best fit still differs from a Drude form for the conductivity. 

After some trial and error we came up with = 0.003 meV and c = 0.2. Our new 
results which include inelastic as well as impurity scattering are shown as the solid triangles 
in Fig. |1]. The agreement with the data is clearly greatly improved over the pure limit, 
i.e.: F"*" = 0, in Eq. (|l]), particularly at low temperatures and for the smaller microwave 
frequencies used in the experiment. The over all fit to the entire data set is quite good but 
certainly not excellent. The theoretical calculations do reproduce well the general trends 
such as the decrease in peak height with increasing microwave frequency and its shift to 
higher energies. The inelastic scattering largely controls this trend and the present analysis 
provides support for the validity of the charge carrier-fluctuation spectrum spectral weight 
Px{^) used here. We stress that the form of Px{^) "was not adjusted in any way to fit the 
microwave data but comes from consideration of the infrared optical scattering rate only. 

One can further examine the quality of the fit by plotting the same data as a function 
of frequency for fixed temperature T. This is done in Fig. |^. In the figure the black solid 
symbols are experiment, gray solid theory with impurities (F"*" = 0.003 meV, c = 0.2) and 
the open symbols the pure case, i.e.: including only the inelastic scattering captured in our 
model spectral density /^x(u;). It is quite clear that some impurity scattering is needed to 
get even reasonably close to the data although a tight fit is never possible. The data do not 
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vary as rapidly at the lowest frequency values as in theory. We point out, however, that 
the fit is very much better than we could achieve using a pure BCS formalism. Inclusion 
of inelastic scattering is essential in any serious attempt to understand this data even at 
reasonably low temperatures. Our use of the generalized Eliashberg equations, given in the 
previous section, can be viewed as a phenomenology with kernel J^x(c<j) determined from 
experimental data. We have previously found that this phenomenology is able to explain 
many of the anomalous superconducting properties observed in the oxides. The present 
study extends the range of agreement to the microwave data although some discrepancies 
do remain. These are not large, however. In Fig. ^ the data for 20 K fit perfectly, for lOK we 
have only disagreement for 1 GHz and at 15 K we have a slight disagreement for 1 and 2 GHz. 
But for these two temperatures Hosseini et a/.0 also have problems with their Drude fits. 
The over all fit seems to be as good as the one of Hosseini et a/.0 and also as in the analysis 
of Berlinsky et a/.0 who conclude that the data does not support a quasiparticle picture. 
Here we find, instead, no serious disagreement of the data with an Eliashberg formulation 
of the d-wave state which includes some impurity scattering described with an intermediate 
value of c, the parameter that spans the interaction strength from unitary (c = 0) to Born 
(c = oo). 

IV. CONCLUSION 

We have made use of a set of charge carrier-boson exchange spectral densities I'^xi^) 
obtained previously from an analysis of infrared optical conductivity data to calculate the 
microwave conductivity at several frequencies as a function of temperature in a generalized 
Eliashberg formalism suitable to describe a ci-wave superconductor. Agreement with recent 
data on pure samples of YBa2Cu306.99 is satisfactory provided a small amount of elastic 
impurity scattering is also included. The impurity potential used is neither in the Born 
(weak) nor unitary (strong) limit. A potential of intermediate strength is indicated. The 
low temperature behavior found in the theory cannot accurately be described by a Drude 
form and does not support the use of a two fluid model with the normal component described 
by a scattering rate constant in frequency although temperature dependent. Instead, the 
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conductivity as a function of energy at fixed T is concave upward reflecting the intrinsic 
frequency dependence of the combined scattering rates. This holds even when inelastic 
scattering is included. The calculations show clearly that, on the hold, the main features 
of the microwave data can be understood within the same generalized Eliashberg formalism 
that has recently been so successful in describing many of the anomalous superconducting 
state propertiesi seen in the oxides. This follows also from the observation by Schachinger 
and Carbotte@ that adding elastic impurity scattering only affects the low energy region 
of the optical properties while the inelastic scattering effects are seen in the energy region 
> Aq. Thus, adding elastic impurity scattering allows a fit of theoretical Eliashberg results 
to match low energy optical properties of a particular sample without violating all earlier 
findings which particularly concentrated on the energy region a; > Aq or on bulk effects. 
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FIG. 1: Microwave conductivity ai{u),t) in 10''il^^m~^ vs. the reduced temperature t = T/Tc 
for the five frequencies measured in experimental work of Hosseini et alE^ namely 17 = 1.14 GHz, 
2.25,13.4, 22.7, and 75.3 GHz (bottom frame). Solid squares are experiment, open triangles clean 
limit and solid triangles inelastic scattering plus impurities characterized by a potential with r"*" = 
0.003 meV and c = 0.2. 

FIG. 2: The self consistent impurity scattering rate 7(c) in the superconducting d-wave state given 



by Eq. (10) for various values of the impurity potential c. It is largest for c = which corresponds 
to the unitary limit and rapidly becomes small as c increases beyond 0.3. The corresponding 
normal scattering rate that we should compare 7(c) with is 7rr+/(l + c^). Here, r+ = 0.15 meV 
and Ao = 24^2 meV. 



FIG. 4: The real part of the optical conductivity cri(uj) in BCS theory with F"*" = 0.15 meV at 
temperature T = lOK for various values of the impurity parameter c, namely Born limit (c 00, 
solid), unitary limit (c = 0, dotted), c = 0.4 (dashed), and c = 1.0 (dash-dotted). The solid gray 
curve is for comparison and gives the normal state. The inset shows the same results on a different 
vertical scale. 



FIG. 5: The microwave conductivity ai{uj,T) as a function of to for three different temperatures. 
The data is the same as shown in Fig. |l]. The open symbols are theory for the pure limit, the 
solid gray symbols theory with some impurity scattering additionally included, and the solid black 
symbols are experiments. The squares are for T = 10 K, the upward triangles for T = 15 K, and 
the down triangles for T = 20K. 



FIG. 3: The imaginary part of the renormalized frequency 0)2(1^) as a function of u for several values 
of the impurity parameter c, namely c = (solid), c = 0.1 (dashed), c = 0.2 (dotted), c = 0.3 
(dash-dotted), and c = 0.4 (dash-double-dotted). Here r+ = 0.15 meV and Aq = 24\/2meV. 
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